Single ion heat engine with maximum efficiency at maximum power 
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We propose an experimental scheme to realize a nano heat engine with a single ion. An Otto cycle 
may be implemented by confining the ion in a linear Paul trap with tapered geometry and coupling it 
to engineered laser reservoirs. The quantum efficiency at maximum power is analytically determined 
in various regimes. Moreover, Monte Carlo simulations of the engine are performed that demonstrate 
its feasibility and its ability to operate at maximum efficiency of 30% under realistic conditions. 
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Miniaturization has lead to the development of increas- 
ingly smaller devices. This ongoing size reduction from 
the macroscalc to the nanoscale is approaching the ul- 
timate limit, given by the atomic nature of matter [1]. 
Prominent macro-devices are heat engines that convert 
thermal energy into mechanical work, and hence motion 
[2]. A fundamental question is whether these machines 
can be scaled down to the single particle level, while re- 
taining the same working principles as, for instance, those 
of a car engine. It is interesting to note in this context 
that biological molecular motors are based on completely 
different mechanisms that exploit the constructive role of 
thermal fluctuations [3, 4]. At the nanoscale, quantum 
properties become important and have thus to be fully 
taken into account. Quantum heat engines have been the 
subject of extensive theoretical studies in the last fifty 
years [5-14]. However, while classical micro heat engines 
have been fabricated, using optomechanical [15], micro- 
electromechanical [16-18], and colloidal systems [19], to 
date no quantum heat engine has been built. 

In this paper, we take a step towards that goal by 
proposing a single ion heat engine using a linear Paul 
trap. Specifically, we present a scheme which has the 
potential to implement a quantum Otto cycle using 
currently available state-of-the-art ion-trap technology. 
Laser-cooled ions in linear Paul traps are quantum sys- 
tems with remarkable properties [20]: they offer an un- 
precedented degree of preparation and control of their 
parameters, permit their cooling to the ground state, and 
allow the coupling to engineered reservoirs [21]. For these 
reasons, they have played a prominent role in the exper- 
imental study of quantum computation and information 
processing applications [22, 23]. They are also invaluable 
tools for the investigation of quantum thermodynamics 
[24]. The quantum Otto cycle for a harmonic oscillator 
is a quantum generalization of the common four-stroke 
car engine and a paradigm for thermodynamic quantum 
devices [25-27]. It consists of two isentropic processes 
during which the frequency of the oscillator (the trap 
frequency) is varied, and of two isochoric processes, that 



in this situation correspond to a change of temperature 
at constant frequency, see Fig. 1(a). In the present pro- 
posal, we simulate the Otto cycle by confining a single ion 
in a novel trap geometry with an asymmetric electrode 
configuration (see Fig. 1(c)) and coupling it alternatingly 
to two engineered laser reservoirs. As for all realistic ma- 
chines, this Otto engine runs in finite time and has thus 
non-zero power [28] . We determine the efficiency at max- 
imum power in the limit of adiabatic and strongly nona- 
diabatic processes, which we express in terms of the nona- 
diabaticity parameter introduced by Husimi [29] . We fur- 
ther present semiclassical Monte Carlo simulations, with 
realistic parameters, that demonstrate the experimental 
feasibility of such a device. The proposal and the sin- 
gle ion trap design idea have several unique advantages: 
First, all of the parameters of the engine, in particular the 
temperatures of the baths, are tunable over a wide range, 
in contrast to existing engines. As a result, maximum ef- 
ficiency can be achieved. Moreover, at low temperatures, 
the engine may operate in the quantum regime, where 
the discreteness of the energy spectrum plays an impor- 
tant role. In addition, the coupling to the laser reservoirs 
can be cither switched on and off externally, or by the in- 
trinsic dynamics of the ion itself. In this latter mode, the 
heat engine runs autonomously [30]. Since trapped ions 
are perfect oscillator models, the results described here 
may in principle be extended to analogous systems, such 
as micro- and nanomechanical oscillators [31-34], offering 
a broad spectrum of potential applications. 

Quantum Otto cycle. We consider a quantum engine 
whose working medium is a single harmonic oscillator 
with time-dependent frequency w t , changing between u\ 
and u>2- The engine is alternatingly coupled to two heat 
baths at inverse temperatures = l/^sTj) (i = 1,2), 
where ks is the Boltzmann constant. The Otto cycle 
consists of four consecutive steps as shown in Fig. 1(a): 
(1) Isentropic compression A(wi,/?i) — > B(u} 2 , fli)- the 
frequency is varied during time n while the system is 
isolated. The evolution is unitary and the von Neumann 
entropy of the oscillator is thus constant. Note that state 
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B is non-thermal even for slow (adiabatic) processes. 

(2) Hot isochore B(w 2 ,/?i) —> C(cj 2 , ^2): the oscillator is 
weakly coupled to a reservoir at inverse temperature /3 2 
at fixed frequency and allowed to relax during time t 2 to 
the thermal state C. This equilibration is much shorter 
than the expansion/compression phases (see below). 

(3) Isentropic expansion C(w 2 ,/3 2 ) — > D(uj\, P2): the fre- 
quency is changed back to its initial value during time 
t 3 . The isolated oscillator evolves unitarily into the non- 
thermal state D at constant entropy. 

(4) Cold isochore -D(wi,,3 2 ) —> A(uj\, fi-i)'. the system 
is weakly coupled to a reservoir at inverse temperature 
Pi > /?2 and quickly relaxes to the initial thermal state 
A during t&. The frequency is again kept constant. 

In order to determine the efficiency of the quantum 
Otto cycle, we need to evaluate work and heat for each of 
the above steps. During stroke (2) and (4), the frequency 
is constant, and thus only heat is exchanged with the 
reservoirs. On the other hand, during stroke (1) and 
(3), the system is isolated and only work is performed by 
modulating the frequency. Since the dynamics is unitary 
in the latter, the Schrodinger equation for the parametric 
oscillator can be solved exactly and its mean energy can 
be obtained analytically. The average quantum energies 
(H) of the oscillator at the four stages of the cycle are 
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FIG. 1. (Color online) (a) Energy-frequency diagram of the 
Otto cycle for the radial degree of freedom of the ion. The 
continuous line represents the ideal process, while the dots 
show the results of the Monte Carlo simulations, (b) The 
pictograms illustrate the four individual strokes of the cycle 
for the radial thermal state, (c) Geometry of the tapered Paul 
trap: the rf-electrodes have an angle of 9 = 20° with the trap 
axis, the length of the trap is 5 mm, the radial distance of the 
ion to the rf-electrodes is ro = 1 mm. 
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where we have introduced the two adiabaticity parame- 
ters Ql and Q 2 [29] . They are equal to one for adiabatic 
(slow) processes and increase with the degree of nonadia- 
baticity. Their explicit expressions for any given modula- 
tion u>t, can be found in Refs. [35, 36]. Equations (la-d) 
reduce to their classical limits when h 0. The mean 
work, denoted by (W\), done during the first stroke is 



and fourth stroke are given by 
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For a heat engine, heat is absorbed from the hot reservoir, 
{Q2) > 0, and flows into the cold reservoir, (Q4) < 0. As 
a result, the two conditions have to be satisfied: 
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whereas the mean heat (Q2) exchanged with the hot 
reservoir during the second stroke reads, 



(Q2) = 



\H)c- 



huj 2 

~2~ 



H) 



coth 



B 



Qi coth 



(3) 



In a similar way, the average work and heat for the third 



The efficiency of this quantum engine, defined as the ratio 
of the total work per cycle and the heat received from the 
hot reservoir, then follows as 
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The above exact expression is valid for arbitrary tem- 
peratures and frequency modulations, and allows for a 
detailed investigation of the performance of the engine. 
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Efficiency at maximum power. Two essential charac- 
teristics of a heat engine are the power output, P = 
-((Wi) + (W 3 ))/(ti + t 2 + t 3 + t 4 ), and the efficiency 
at maximum power [28]. Both can be evaluated ana- 
lytically for the quantum Otto cycle with the help of 
Eq. (7). We shall separately consider the case of an adi- 
abatic process, Q\ 2 — T an d the case of a sudden switch 
of the frequencies for which Q* 2 = (ujf + uj 2 )/(2ojiuj 2 ). 
Let us begin with the high temperature regime fithujj <C 
= 1,2). The total work produced by the heat 
engine for a quasistatic frequency modulation is given by 
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Assuming that the initial frequency lui is fixed and by 
optimizing with respect to the second frequency W2, we 
find that the power is maximum when uj 2 /^>i — \J Pi/ fi 2 - 
As a consequence, the efficiency at maximum power is 
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which corresponds to the Curzon-Ahlborn efficiency [37]. 
Conversely, for a sudden switch, the total work is 
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By optimizing again with respect to lu 2 , we find the power 
to be maximized when the condition ui 2 /uji — \j fli is 
satisfied. The corresponding efficiency reads [26], 
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Equations (9) and (11) show that maximum efficiency 
(of either 1 or 1/2) can be attained when j3i — > oo. In 
this low-temperature limit, quantum effects are crucial. 
Repeating the above optimization analysis in the regime 
f>ihu>i 3> 1, we find the efficiency at maximum power 
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for a sudden frequency switch when uj 2 = \j2iJ\lhji 2 . 
The above expressions, in which the classical thermal en- 
ergy ksTi is replaced by the ground state energy hu!i/2 of 
the oscillator, are the quantum extensions of the Curzon- 
Ahlborn and Rczck-Kosloff results (9) and (11). 

Proposed realization in a Paul trap. Such single ion 
heat engine is composed of one trapped ion in a modified 
linear Paul trap, as sketched in Fig. 1(c). The trapped 
ion is initially prepared in a thermal state at low temper- 
ature by laser cooling to the Doppler limit in all spatial 
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FIG. 2. (Color online) Efficiency at maximum power of the 
Otto engine as a function of the temperature ratio, (a) shows 
the Carnot efficiency with zero power, (b) corresponds to the 
theoretical adiabatic process, Eq. (9), while (c) corresponds 
to the sudden frequency switch, Eq. (11). The black points 
denote the results of the numerical simulations and demon- 
strate that the engine runs at maximum efficiency at maxi- 
mum power. The horizontal bars indicate estimated upper 
bounds at CJ2 = 1.5 wi. 



directions. The engine is driven by alternatingly coupling 
the ion to two reservoirs that heat and cool the thermal 
state of the ion in the radial direction through scattering 
forces. These two baths are realized by blue and red de- 
tuned laser beams on a cycling transition of the trapped 
ion, irradiated in the radial plane (x,y). In a first step, 
the coupling to the heat reservoirs is switched on and off 
externally. The geometry of a trap design with tapered 
radio-frequency (rf) electrodes leads to a pseudopotential 
of the form [39] , 
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where is the angle between the electrodes and the 
trap axis z, and ro the radial distance of the ion to 
the electrodes, as shown in Fig. 1(c). This potential re- 
sults in radial trap frequencies that depend on the ax- 
ial position z, and in an axial force that depends on 
the radial displacement. The coupling between har- 
monic axial and radial modes is of the generic form, 
H = T.teix^y^z} huj 0t (1/2 + aj ai)~Cz- {uj% x x 2 + wg y y 2 ) , 
valid for small z, where C = 2 m tan 0/ro denotes the cou- 
pling constant between the oscillator modes, and woi are 
the trap frequencies at the center of the trap. A change 
in temperature T of the radial thermal state of the ion, 
and thus of the width of its spatial distribution, leads to a 
modification in the axial component of the repelling force 
which changes the point of equilibrium zq(T). Heating 
and cooling the radial state hence moves the ion back 
and forth along the trap axis, as sketched in Fig. 1(b), 
resulting in the closed Otto cycle shown in Fig. 1(a). 
This thermally induced axial movement corresponds to 



the mechanically usable movement of a piston of a classi- 
cal engine, while the radial mode corresponds to the gas 
in the cylinder. 

The energy gained by running the engine cyclically in 
the radial direction can be stored in the axial degree of 
freedom. Indeed, the axial translation on the nanometer 
scale can be enhanced by coupling the ion to the laser 
reservoirs resonantly at the axial trap frequency. In such 
a way, the cyclic temperature change of the radial ther- 
mal state is directly converted into an increasing coherent 
axial oscillation. In principle, the excited axial oscillation 
is only limited by the trap geometry. The cycle time of 
the engine is given by the axial oscillation period, and 
is of the order of 10 /is, while the time needed to change 
the temperature is below 1 /is. In order to reach a steady 
state of the heat engine, a red-detuned low intensity dis- 
sipation laser is applied in axial direction that damps 
the coherent movement depending on the velocity of the 
ion. The stored energy in the axial motion may be trans- 
ferred to other oscillator systems, e.g. separately trapped 
ions [38] or nanomechanical oscillators [40], and thus ex- 
tracted from the engine as usable work. When driven as 
a heat pump, cooling of such systems should be possible. 

The single ion engine offers complete control of all pa- 
rameters over a wide range. The temperature of the two 
heat reservoirs, as well as the applied dissipation, can 
be adjusted by tuning the laser frequencies and intensi- 
ties. On the other hand, the oscillation amplitude and 
frequency of the ion can be changed by virtue of the trap 
parameters. This unique flexibility of the device can be 
exploited to satisfy the conditions for maximum efficiency 
at maximum power derived previously. 

Monte Carlo simulations. We have performed exten- 
sive semiclassical simulations of the engine using Monte- 
Carlo and partitioned Runge-Kutta integrators, as de- 
scribed in [41]. We reach a dynamic confinement of the 
ion in the trap volume through an oscillating parabolic 
and tapered saddle potential of the form 
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with a trap drive frequency w r f = 60 MHz, and the trap 
geometry described in Fig. 1(c). The thermal state of 
the ion is generated by a Boltzmann distributed ensem- 
ble over several thousand classical trajectories [42]: the 
initial parameters for each ion are choosen randomly, and 
the desired thermal probability distribution is reached 
through the ion-light interaction and the corresponding 
stochastic spontaneous emission of photons [43]. The 
switching of the detuned heating and cooling lasers is 
adjusted to the axial trap frequency uj z . Each laser is 
coupled to the ion for 20% of an axial trap period. The 
ensemble of driven oscillators is thus excited coherently 
in axial direction such that heating/cooling takes place at 
the turning points of the trajectory, independently of the 
random initial conditions. In phase space, the ensemble 




FIG. 3. (Color online) Motional-state phase space distribu- 
tion during an engine cycle ABCD: (a) radial mode changing 
from thermal to nonthermal states, (b) In axial direction, a 
coherent oscillation with typical displacement of a — 10 3 wave 
packets displays the functional performance of the engine. 



average of the axial mode performs a phase synchronized 
oscillation (see Fig. 3), while the radial mode changes 
between thermal and nonthermal states. 

Radial temperatures in the range of 20 to 200 mK were 
achieved, corresponding to 0.1 < 02/01 < 1 and respec- 
tive radial phonons numbers of about 400 and 4000. For 
a realistic maximal axial amplitude of about 1 mm, the 
relative variation of the radial frequency at 6.0 MHz is 
about 50%. By properly adjusting the parameters to 
satisfy the optimality condition, our simulations show 
that this Otto engine has the ability to run at maxi- 
mum efficiency at maximum power in the interval 0.5 < 
02 1 01 < 1 (see Fig. 2). The maximum efficiency, deter- 
mined through Eq. (7) with values taken from simulated 
cycles as shown in Fig. 1(a), is about 30%. This maxi- 
mum efficiency is significantly larger than those obtained 
to date [18, 19]. The expected power of the engine is of 
the order of P ~ 10 _20 J/s. The case of a sudden switch 
can be realized by exchanging the values of radial and ax- 
ial frequencies. Contrary to the adiabatic limit, the op- 
timality condition can be achieved down to 02 / 01 — 0.2. 
However, the maximum efficiency reduces to 23%. 

While so far the laser reservoirs have, for simplicity, 
been switched on and off externally, another feature of 
this Otto engine is the ability to operate in a completely 
self-driven manner. To this end, the foci of the heating 
and cooling lasers can be separated spatially on the trap 
axis by e.g. 200 fim, so that the ion is coupled to the 
heat baths at the turning points of its axial trajectory. 
No active switching is required and the axial motion of 
the ion is self-amplifying. The ion needs to be driven 
only in the initial phase of the axial motion to reach a 
threshold amplitude. We note that at a radial frequency 
of 6.0 MHz, the thermal energy of the oscillator, Eq. (la), 
starts to deviate appreciably from its classical value be- 
low 200 /iK, which could be reached with 40 Ca-ions, if we 
assume a two-level approximation and the Doppler cool- 
ing limit Td — HT / {2ks), where T denotes the linewidth 
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of the dipolc transition. A single ion engine has thus 
the potential to enter the quantum regime and become 
a tool to study effects of quantum coherence and corre- 
lations on the efficiency. Application of optimal control 
techniques [44] would further allow for nonclassical bath 
engineering. The investigation of heating and cooling 
on the simple and fundamental single ion mode interac- 
tions may serve for prototyping coolers also in systems, 
which share similar properties. For the specific example 
of micromechanical oscillators, mode coupling has been 
described and realized in several experiments [33, 34]. 

Conclusion. We have put forward a realistic proposal 
for a tunable nanoengine based on a single ion in a ta- 
pered linear Paul trap coupled to engineered laser reser- 
voirs. The operation in the Otto cycle would result in 
coherent ion motion. Combining analytical and numeri- 
cal analysis, we have studied the performance of the en- 
gine and showed that it can achieve maximum efficiency 
at maximum power in a wide range of temperatures. We 
expect a stimulating impact on the development on nano- 
and micromechanical oscillators which are of fundamen- 
tal importance for future sensor technologies. 
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